Sensor coordinate calibration in an ultrasound system

ABSTRACT

A medical imaging apparatus includes: an ultrasound probe to transmit and receive an ultrasound signal to obtain an ultrasound image of an object; a sensor coupled to the ultrasound probe to provide a position information related to an ultrasound imaging position in the object; and a processor to extract a certain anatomical feature from the ultrasound image and the certain anatomical feature from a medical image, perform a first registration between the medical image and the ultrasound image by matching the certain anatomical feature in the ultrasound image the certain anatomical feature in the medical image, perform a second registration between the sensor and the medical image based on the position information and the first registration, and control a display to display a portion of the medical image corresponding to the ultrasound imaging position based on at least one of the first registration and the second registration.

CROSS-REFERENCE TO RELATED APPLICATIONS

The present application is a Continuation of U.S. patent applicationSer. No. 15/943,061, filed Apr. 2, 2018, which is a Continuation of U.S.application Ser. No. 14/752,165, filed Jun. 26, 2015, now U.S. Pat. No.9,955,951 issued May 1, 2018, which is a Continuation of U.S.application Ser. No. 12/846,545, filed Jul. 29, 2010, now U.S. Pat. No.9,082,178 issued Jul. 14, 2015, which claims priority from Korean PatentApplication No. 10-2009-0070994, filed Jul. 31, 2009. The disclosure ofeach of these applications is incorporated herein by reference in itsentirety.

TECHNICAL FIELD

The present invention generally relates to ultrasound systems, and moreparticularly to an ultrasound system and method for performing sensorcoordinate calibration through image-based registration between athree-dimensional ultrasound image and a computed tomography (CT) image.

BACKGROUND

The ultrasound system has become an important and popular diagnostictool due to its non-invasive and non-destructive nature. Modernhigh-performance ultrasound imaging diagnostic systems and techniquesare commonly used to produce two- or three-dimensional images ofinternal features of patients (target objects).

However, the ultrasound system suffers from inherent shortcomings of anultrasound image such as a low signal-to-noise ratio and a limited fieldof view. Thus, the image registration of a CT (or MR) image onto theultrasound image has been introduced in order to compensate fordeficiencies of the ultrasound image. A sensor has been used to performthe image registration of a CT (or MR) image onto the ultrasound image.Researches have been introduced to calibrate the sensor to matchcoordinates of the CT image and coordinates of the sensor.

Conventionally, after outer markers are attached on a surface of atarget object, a CT image and an ultrasound image for the target objectwith the markers are acquired. Thereafter, the calibration is carriedout by using a relationship between coordinates of the markers in the CTand ultrasound image. That is, the outer markers should be attached tothe surface of the target objects before obtaining the CT image and theultrasound image and be maintained in the same position until completingthe acquisition of the ultrasound image. Moreover, a sensor must sensethe positions of the respective outer markers.

Further, the registration between the coordinate of the CT image and thecoordinate of the sensor has been performed by manually inputting innermarkers on the CT image. Thus, a user of the ultrasound system had toinput the inner markers, which causes the registration between thecoordinate of the CT image and the coordinate of the sensor to be wrong.

SUMMARY

An embodiment for forming a plurality of three-dimensional ultrasoundimages is disclosed herein. In one embodiment, by way of non-limitingexample, an ultrasound system, includes: an ultrasound image formingunit including a ultrasound probe and being configured to form athree-dimensional ultrasound image of a target object; a sensor coupledto the ultrasound probe; a memory configured to store athree-dimensional computed tomography (CT) image of the target objectand position information on a position between the three-dimensionalultrasound image and the sensor; and a processor configured to performimage registration between the three-dimensional CT image and thethree-dimensional ultrasound image to thereby form a firsttransformation function for transforming a position of the sensor to acorresponding position on the three-dimensional CT image and performcalibration of the sensor by applying the position information to thefirst transformation function.

In another embodiment, a method of performing a calibration of a sensor,includes: a) obtaining a three-dimensional ultrasound image of a targetobject obtained by the ultrasound system and a three-dimensional CTimage; b) calculating a position information on a position between thethree-dimensional ultrasound image and the sensor; c) performingregistration between the three-dimensional ultrasound image and thethree-dimensional CT image to obtain a first transformation function fortransforming a position of the sensor to a corresponding position on thethree-dimensional CT image; and d) performing calibration of the sensorby applying the position information to the first transformationfunction.

The Summary is provided to introduce a selection of concepts in asimplified form that are further described below in the DetailedDescription. This Summary is not intended to identify key or essentialfeatures of the claimed subject matter, nor is it intended to be used indetermining the scope of the claimed subject matter.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a block diagram showing an illustrative embodiment of anultrasound system.

FIG. 2 is a block diagram showing an illustrative embodiment of anultrasound image forming unit.

FIG. 3 is an illustrative embodiment of an ultrasound probe.

FIG. 4 is a block diagram showing an illustrative embodiment of aprocessor.

FIG. 5 is a schematic diagram showing an example of eigenvalues in theHessian matrix.

DETAILED DESCRIPTION

A detailed description may be provided with reference to theaccompanying drawings. One of ordinary skill in the art may realize thatthe following description is illustrative only and is not in any waylimiting. Other embodiments of the present invention may readily suggestthemselves to such skilled persons having the benefit of thisdisclosure.

FIG. 1 is a block diagram showing an illustrative embodiment of anultrasound system. The ultrasound system 100 may include an ultrasoundimage forming unit 110, a sensor 120, a memory 130, a processor 140 anda display unit 150.

The ultrasound image forming unit 110 may be configured to transmitultrasound signals to a target object (not shown) and receive ultrasoundecho signals reflected from the target object. The ultrasound imageforming unit 110 may be further configured to form a three-dimensionalultrasound image of the target object based on the received ultrasoundecho signals.

FIG. 2 is a block diagram showing an illustrative embodiment of anultrasound image forming unit 110. The ultrasound image forming unit 110may include a transmit (Tx) signal generating section 111, an ultrasoundprobe 112 including a plurality of transducer elements (not shown), abeam former 113, an ultrasound data forming section 114 and an imageforming section 115.

The Tx signal generating section 111 may generate Tx signals accordingto an image mode set in the ultrasound system 100. The image mode mayinclude a brightness (B) mode, a Doppler (D) mode, a color flow mode,etc. In one exemplary embodiment, the B mode may be set in theultrasound system 100 to obtain a B mode ultrasound image.

The ultrasound probe 112 may receive the Tx signals from the Tx signalgenerating section 111 and generate ultrasound signals, which may travelinto the target object. The ultrasound probe 112 may further receiveultrasound echo signals reflected from the target object and convertthem into electrical receive signals. In such a case, the electricalreceive signals may be analog signals. The ultrasound probe 112 may be athree-dimensional probe, a two-dimensional probe, a one-dimensionalprobe or the like.

FIG. 3 is an illustrative embodiment of an ultrasound probe 112. Atleast one transducer element (not shown) of the ultrasound probe 112generates an image plane IP, which is used to scan a region of interestROI. The image plane IP may be one of slice planes of thethree-dimensional ultrasound image. The sensor 120 is attached to thehousing of the ultrasound probe 112 to determine the position andorientation of the image plane IP. The ultrasound system 100 coupledwith the ultrasound probe 112 via the probe cable 105 can use the datagenerated by the sensor 120 to determine the position and orientation ofthe sensor 120 and/or the image plane IP, as described below.

In this preferred embodiment, the sensor 120 is a magnetic sensor thatmonitors the free-hand movement of the ultrasound probe 112 in sixdegrees of freedom with respect to a transducer element 170. As shown inFIG. 3, the sensor 120 and the transducer element 170 each define anorigin (122, 172, respectively) defined by three orthogonal axes (X′,Y′, Z′ and X″, Y″, Z″, respectively). The sensor 120 monitors thetranslation of the origin 122 with respect to the origin 172 of thetransducer element 170 to determine position and monitor the rotation ofthe X′, Y′, Z′ axes with respect to the X″, Y″, Z″ axes of thetransducer element 170 to determine orientation.

The position and orientation of the sensor 120 can be used to determinethe position and orientation of the image plane IP. As shown in FIG. 3,the image plane IP defines an origin OR defined by three orthogonal axesX, Y, Z, which are preferably aligned with the origin of a centeracoustic line generated by the ultrasound probe 112. The position of theorigin 122 and the orientation of axes X′, Y′, Z′ of the sensor 120 maynot precisely coincide with the position of the origin OR and theorientation of the axes X, Y, Z of the image plane IP. For example, inFIG. 3, the origin OR of the image plane IP is offset from the origin122 of the sensor 120 by a distance z₀ along the Z-direction and adistance of y₀ along the Y-direction. In FIG. 3, there is no offsetalong the X-direction nor is there a rotational offset in theorientation of the axes. Accordingly, the position and orientation ofthe sensor 120 do not directly describe the position and orientation ofthe image plane IP.

To determine the position an orientation of the image plane IP from theposition and orientation of the sensor 120, sensor calibration data isused to transform the position and orientation of the sensor 120 to theposition and orientation of the image plane IP. For simplicity, the term“position and orientation” is used to broadly refer to position and/ororientation. Accordingly, if the sensor 120 has the same orientation asthe image plane IP, then the position and orientation calibration datamay not contain any orientation calibration data. Similarly, as shown inFIG. 3, the sensor 120 may not have a positional offset with respect toone or more axes of the image plane IP.

There are a number of ways of defining the image plane/sensor offset.One method of calibrating at least some types of sensors use threeorthogonal linear dimension offsets in X, Y, Z and three rotation anglesabout each of these axes. Other methods include using a positiontransformation matrix or quaternions, which are described in the usermanual for the mini Bird™ and the Flock of Bird™ systems by AscensionTechnology Corp.

As described above, the ultrasound probes with position and orientationsensors are typically used only with ultrasound systems that contain thecalibration data for the probe/sensor pair. Conventionally, theprobe/sensor pair is calibrated, and the calibration data is stored inthe ultrasound system 100, which will be used in conjunction with theprobe/sensor pair. If the probe/sensor pair is to be used with adifferent ultrasound system, then the probe/sensor pair typically needsto be re-calibrated on that different ultrasound system. Sincesonographers are often unable or unwilling to perform probe/sensor paircalibration, probe/sensor pairs are often used only with the ultrasoundsystem for which the probe/sensor pair was initially calibrated.

Referring back to FIG. 2, the beam former 113 may convert the electricalreceive signals outputted from the ultrasound probe 112 into digitalsignals. The beam former 113 may further apply delays to the digitalsignals in consideration of the distances between the transducerelements and focal points to thereby output receive-focused signals.

The ultrasound data forming section 114 may form a plurality ofultrasound data by using the receive-focused signals. In one embodiment,the plurality of ultrasound data may be radio frequency (RF) data or IQdata. The image forming section 115 may form the three-dimensionalultrasound image of the target object based on the ultrasound data.

Referring back to FIG. 1, the sensor 120 may be mounted on one side ofthe ultrasound probe 112. In one embodiment, by way of non-limitingexamples, the sensor 120 may be built in the ultrasound probe 112 to beaway from the plurality of transducer elements (not shown) by apredetermined distance. Alternatively, the sensor 120 may be externallymounted on the ultrasound probe 112 to be away from the plurality oftransducer elements. The sensor 120 may include three-dimensionalsensor, which can detect a three-dimensional position and an angle ofthe ultrasound probe 112.

The memory 130 may store a three-dimensional CT image of the targetobject. In one embodiment, by way of non-limiting examples, thethree-dimensional CT image may be a three-dimensional CT image of aliver in which a diaphragm and a blood vessel are extracted. The memory130 may store information on a position between the three-dimensionalultrasound image and the sensor 120 (hereinafter, referred to as“position information”). The position information may includeinformation on a distance between the transducer elements (not shown)and the sensor 120. In one embodiment, by way of non-limiting examples,the memory 120 may include at least one of a random access memory (RAM),a hard disk drive or the like.

The processor 140 may be configured to perform registration between thethree-dimensional CT image and the three-dimensional ultrasound image,thereby forming a transformation function (T_(probe)) for representingthe ultrasound probe 112 on the three-dimensional CT image. Furthermore,the processor 140 may perform calibration of the sensor 120 to matchcoordinates of the three-dimensional CT image (not shown) andcoordinates of the sensor 120 based on the position information and thetransformation function.

FIG. 4 is a block diagram showing an illustrative embodiment of theprocessor 140. The processor 140 may include a diaphragm extractingsection 141, a vessel extracting section 142, a diaphragm refiningsection 143, a registration section 144, a calibration section 145 andan image processing section 146.

The diaphragm extracting section 141 may be configured to extract adiaphragm from the three-dimensional ultrasound image formed in theultrasound image forming unit 110. In one embodiment, the diaphragmextracting section 141 may perform a Hessian matrix based flatness testupon the three-dimensional ultrasound image to extract the diaphragm.The diaphragm may be considered as a curved surface in thethree-dimensional ultrasound image. Thus, regions in which a voxelintensity change in a normal direction at a surface is greater than avoxel intensity change in a horizontal direction at the surface may beextracted as the diaphragm. FIG. 5 is a schematic diagram showing anexample of eigenvalues λ1, λ2 and λ3 in the Hessian matrix.

Hereinafter, an operation of the diaphragm extracting section 141 willbe described in detail. The diaphragm extracting section 141 may beconfigured to select voxels having a relatively high flatness value. Theflatness μ(v) may be defined as the following equation (1).

μ(v)=ϕ₁(v)ϕ₂(v)ϕ₃(v)/ϕ₃ _(max) (v)   (1)

wherein ϕ₁(v), ϕ₂(v) and ϕ₃(v) in the equation (1) may be represented asthe following equation (2).

$\begin{matrix}{{{\varphi_{1}(v)} = ( {1 - \frac{\lambda_{1}(v)}{\lambda_{3}(v)}} )^{2}},{{\varphi_{2}(v)} = ( {1 - \frac{\lambda_{2}(v)}{\lambda_{3}(v)}} )^{2}},{{\varphi_{3}(v)} = {\sum\limits_{i}{\lambda_{i}(v)}^{2}}}} & (2)\end{matrix}$

wherein λ₁, λ₂ and λ₃ denote eigenvalues of the Hessian matrix at voxelv. The flatness μ(v) may be normalized to have values of ranging 0-1. Aflatness map may be formed based on the flatness values obtained fromall of the voxels according to the equations (1) and (2). Thereafter,the voxels having a relatively high flatness value are selected. In oneembodiment, the voxels having the flatness over 0.1 may be selected.

The diaphragm extracting section 141 may be further configured toperform the morphological opening upon the selected voxels to removesmall clutters therefrom. The morphological opening may be carried outby sequentially performing erosion and dilation. That is, apredetermined number of the voxels are removed in the edges of the areain which the voxels exist, and thus, the area becomes contracted(erosion). In this manner, it becomes possible to remove small clutters.Thereafter, the edges of the area are expanded by the predeterminednumber of the voxels (dilation). These erosion and dilation may beperformed by one or more voxels.

The diaphragm is the largest surface in the three-dimensional ultrasoundimage. The largest surface may be selected among candidate surfacesobtained by the intensity-based connected component analysis (CCA) forthe voxels and the selected surface may be regarded as the diaphragm inthe three-dimensional ultrasound image. The voxel-based CCA is one ofthe methods of grouping regions in which voxel values exist. Forexample, the number of voxels connected to each of the voxels through aconnectivity test by referring to values of voxels neighboring thecorresponding voxel (e.g., 26 voxels) may be computed. The voxels, ofwhich connected voxels are greater than the predetermined number, areselected as candidate groups. Since the diaphragm is the widest curvedsurface in the ROI, the candidate group having the most connected voxelsmay be selected as the diaphragm. The surface of the diaphragm may besmoothened.

The vessel extracting section 142 may be configured to perform vesselextraction upon the three-dimensional ultrasound image. The vesselextracting section 142 may be configured to perform the vesselextraction through ROI masking, vessel segmentation and classification.

To avoid mis-extraction of the vessels due to mirroring artifacts, theROI masking may be applied to the three-dimensional ultrasound image bymodeling the diaphragm as a polynomial curved surface. For example, theROI masking may be used to model the diaphragm as the polynomial curvedsurface by using the least means square. However, if all of the lowerportions of the modeled polynomial curved surface are eliminated, thenmeaningful vessel information may be lost at a portion of regions due toan error of the polynomial curved surface. To avoid losing the vesselinformation, the lower portion of the modeled polynomial curved surfacemay be eliminated with a marginal distance. For example, the marginaldistance may be set to about 10 voxels at a lower portion of the ROImask.

Subsequently, the vessel extracting section 142 may be furtherconfigured to segment a vessel region and a non-vessel region. Toexclude non-vessel high intensity regions such as the diaphragm andvessel walls, a low intensity bound value having a less reference boundvalue in the ROI masked three-dimensional ultrasound image may be set asa reference bound value. Thereafter, voxels with a higher intensityvalue than the reference bound value may be removed. The remainingregions may be binarized by using an adaptive threshold value. Then, thebinarized segments may be labeled as vessel candidates.

Next, the vessel extracting section 142 may be further configured toremove non-vessel-type clutters from the binarization image to form realvessel regions from the vessel candidates. In one embodiment, the vesselclassification may include a size test, which evaluates the goodness offit to a cylindrical tube, for filtering out tiny background clutters, astructure-based vessel test for removing non-vessel type clutters, i.e.,an initial vessel test, a gradient magnitude analysis, and a finalvessel test for precisely removing the clutters. Although some cluttersare not removed through the structure-based vessel test, an initialthreshold may be marginally set so that all vessels may be included. Forexample, a threshold value of the initial vessel test may be set to 0.6.At the final vessel test, clutters, which may be formed by small shadingartifacts having low gradient magnitudes, may be precisely removed byconsidering variation of voxel values, i.e., gradient magnitudes, tothereby extract vessel data. In one embodiment, a threshold of the finalvessel test may be set to 0.4.

The diaphragm refining section 143 may be configured to refine thediaphragm region by removing the clutters with the extracted vesselregions. The clutters are mainly placed near the vessel walls.Especially, the vessel walls of inferior vena cava (IVC) are more likelyto be connected to the diaphragm and cause clutters. These clutters maydegrade the accuracy of the feature based registration, and thus, it maybe necessary to refine the diaphragm region. To refine the diaphragm,the vessel regions are extracted according to the vessel extractionmentioned above, the extracted vessel regions may be dilated, and thenthe dilated vessel regions may be subtracted from the initiallyextracted diaphragm region to estimate vessel walls. The estimatedvessel walls may be removed from the diaphragm region. Finally, thediaphragm region may be extracted by applying CCA and the size test.

The registration section 144 may be configured to perform the imageregistration between the three-dimensional ultrasound and CT image. Theregistration section 144 may extract sample points from the vesselregions and the diaphragm region, respectively, among the featuresextracted from the three-dimensional ultrasound image. Also, after thevessel regions and the diaphragm region are extracted from the CT image,the registration section 144 may extract sample points from the vesseland the diaphragm region, respectively. The image registration betweenthe three-dimensional ultrasound and CT image may be performed based onthe extracted sample points to thereby form the transformation function(T_(probe)) between the three-dimensional CT image and thethree-dimensional ultrasound image. The transformation function(T_(probe)) may be given by a matrix and used to transform a position ofthe ultrasound probe 112 to a corresponding position on thethree-dimensional CT image.

The calibration section 145 may perform the calibration of the sensor120 based on the transformation matrix (T_(probe)) from the registrationsection 144 and the position information stored in the memory 130. Moreparticularly, the calibration section 145 may form a transformationmatrix (T_(sensor)) between the sensor 120 and the three-dimensionalultrasound image, i.e., a transformation matrix representing a positionof the sensor 120 with respect to the three-dimensional ultrasoundimage. The transformation matrix (T_(sensor)) may be given by a matrix.The transformation matrix (T_(sensor)) may be defined as the followingequation (3).

$\begin{matrix}{T_{sensor} = {\begin{matrix}{r\; 11} & {r\; 12} & {r\; 13} & x \\{r\; 21} & {r\; 22} & {r\; 23} & y \\{r\; 31} & {r\; 32} & {r\; 33} & z \\0 & 0 & 0 & 1\end{matrix}}} & (3)\end{matrix}$r11=cos θy*cos θz+sin θx*sin θy*sin θz

r12=sin θz*cos θy−sin θx*sin θy*cos θz

r13=cos θx*sin θy, r21=sin θz*cos θx

r22=cos θz*sin θx, r23=sin θx

r31=sin θz*sin θx*cos θy−cos θz*sin θy

r32=−cos θz*sin θx*cos θy−sin θz*sin θy

r33=cos θx*cos θy

wherein, x denotes coordinate of a lateral direction of the sensor 120,y denotes coordinate of an elevation direction of the sensor 120, zdenotes an axial direction of the sensor 120, θ_(x) denotes an angle ofthe sensor 120 from the x-axis, θ_(y) denotes an angle of the sensor 120from the y-axis, and θ_(z) denotes an angle of the sensor 120 from thez-axis. The elevation direction may be a swing direction of thetransducer elements, the axial direction may be a scan line directionfrom the transducer elements and the lateral direction may be alongitudinal direction of the transducer elements.

The calibration section 145 may perform the calibration based on thetransformation matrix (T_(probe)) and the transformation matrix(T_(sensor)). The calibration section 145 may form a transformationmatrix (T) representing the position of the sensor 120 on thethree-dimensional CT image. In one embodiment, the calibration section145 may form the transformation matrix (T) through matrix multiplicationof the transformation matrix (T_(probe)) and the transformation matrix(T_(sensor)).

The image processing section 146 may apply the transformation matrix (T)to the three-dimensional CT image to thereby form a two-dimensional CTimage according to a two-dimensional ultrasound image.

Referring back to FIG. 1, the display unit 150 may display thetwo-dimensional CT image, which is provided from the processor 140.Furthermore, the display unit 150 may display the three-dimensionalultrasound image and the three-dimensional CT image.

Any reference in this specification to “one embodiment,” “anembodiment,” “example embodiment,” “illustrative embodiment,” etc. meansthat a particular feature, structure or characteristic described inconnection with the embodiment is included in at least one embodiment ofthe present invention. The appearances of such phrases in various placesin the specification are not necessarily all referring to the sameembodiment. Further, when a particular feature, structure orcharacteristic is described in connection with any embodiment, it issubmitted that it is within the purview of one skilled in the art toaffect such feature, structure or characteristic in connection withother ones of the embodiments.

Although embodiments have been described with reference to a number ofillustrative embodiments thereof, it should be understood that numerousother modifications and embodiments can be devised by those skilled inthe art that will fall within the spirit and scope of the principles ofthis disclosure. More particularly, numerous variations andmodifications are possible in the component parts and/or arrangements ofthe subject combination arrangement within the scope of the disclosure,the drawings and the appended claims. In addition to variations andmodifications in the component parts and/or arrangements, alternativeuses will also be apparent to those skilled in the art.

What is claimed is:
 1. A medical imaging apparatus comprising: anultrasound probe configured to transmit and receive an ultrasound signalto obtain an ultrasound image of an object; a sensor coupled to theultrasound probe to provide a position information related to anultrasound imaging position in the object; and a processor to extract acertain anatomical feature from the ultrasound image and the certainanatomical feature from a medical image, perform a first registrationbetween the medical image and the ultrasound image based on the certainanatomical feature, perform a second registration between the sensor andthe medical image based on the position information and the firstregistration, and control a display to display a portion of theultrasound image corresponding to the ultrasound imaging position and aportion of the medical image corresponding to the ultrasound imagingposition based on at least one of the first registration and the secondregistration, wherein the processor performs the first registration bymatching the certain anatomical feature in the ultrasound image thecertain anatomical feature in the medical image.
 2. The medical imagingapparatus of claim 1, wherein the processor is further configured toobtain a first sample point of the certain anatomical feature extractedfrom the medical image of the object, obtain a second sample point ofthe certain anatomical feature extracted from the ultrasound image ofthe object, and perform the first registration between the medical imageand the ultrasound image based on the first sample point and the secondsample point.
 3. The medical imaging apparatus of claim 1, wherein theultrasound image is a three-dimensional (3D) ultrasound image, themedical image is a volumetric medical image which includestwo-dimensional (2D) medical images obtained over a region of theobject, the portion of the medical image includes one of the 2D medicalimages, and the processor is further configured to control the displayto display that one of the 2D medical images which corresponds to theultrasound imaging position in the 3D ultrasound image, as indicated bythe sensor.
 4. The medical imaging apparatus of claim 1, wherein theprocessor is further configured to obtain a first registrationinformation which provides a relationship between a coordinate system ofthe medical image and a coordinate system of the ultrasound image basedon the first registration, and calculate a first transform between thecoordinate system of the medical image and the coordinate system of theultrasound image, based on the first registration information, and thefirst transform is applied to perform a conversion between thecoordinate system of the medical image and the coordinate system of theultrasound image.
 5. The medical imaging apparatus of claim 4, whereinthe processor is further configured to perform a third registrationbetween the ultrasound image and the sensor, based on the positioninformation, calculate a second transform which defines a relationshipbetween a coordinate system of the sensor and the coordinate system ofthe ultrasound image, and perform the second registration between themedical image and the sensor based on a combination function of thefirst transform and the second transform.
 6. The medical imagingapparatus of claim 5, wherein the combination function is a product ofthe first transform and the second transform.
 7. The medical imagingapparatus of claim 1, wherein the ultrasound probe comprises one of athree-dimensional ultrasound probe, a two-dimensional ultrasound probe,and a one-dimensional ultrasound probe.
 8. The medical imaging apparatusof claim 1, wherein the sensor comprises a magnetic sensor to sense aposition and an angle of a transducer mounted to the ultrasound probe.9. The medical imaging apparatus of claim 1, wherein the object includesa liver.
 10. The medical imaging apparatus of claim 1, wherein theultrasound image and the medical image are three-dimensional images. 11.The medical imaging apparatus of claim 10, wherein the certainanatomical feature is included in a surface area imaged inside theobject, and the processor is further configured to extract the surfacearea from the medical image of the object, extract a portion of acorresponding surface area from the ultrasound image of the object, andperform the first registration between the medical image and theultrasound image by matching one or more points of the surface areaextracted from the medical image with one or more points of the portionof the corresponding surface area extracted from the ultrasound image.12. The medical imaging apparatus of claim 10, wherein the certainanatomical feature is included in a vessel, and the processor is furtherconfigured to extract the vessel from the medical image of the object,extract a portion of a corresponding vessel from the ultrasound image ofthe object, and perform the first registration between the medical imageand the ultrasound image by matching one or more points of the vesselextracted from the medical image with one or more points of the portionof the corresponding vessel extracted from the ultrasound image.
 13. Themedical imaging apparatus of claim 12, wherein the processor is furtherconfigured to extract the vessel from the medical image of the object bysegmenting the vessel from a region of the object containing a vesselportion and a non-vessel portion in the medical image, and extract theportion of the corresponding vessel from the ultrasound image of theobject by segmenting the portion of the corresponding vessel from theregion of the object containing the vessel portion and the non-vesselportion in the ultrasound image.
 14. The medical imaging apparatus ofclaim 1, wherein the processor is further configured to extract thecertain anatomical feature from the medical image which has beenacquired by a medical imaging modality which captures imaging data viaan imaging technique different from that of an ultrasound apparatus. 15.A medical imaging method comprising: transmitting and receiving anultrasound signal to obtain an ultrasound image of an object; providing,by a sensor of an ultrasound apparatus, a position information relatedto an ultrasound imaging position in the object; extracting a certainanatomical feature from the ultrasound image and the certain anatomicalfeature from a medical image; performing a first registration betweenthe medical image and the ultrasound image, based on the certainanatomical feature; performing a second registration between the sensorand the medical image, based on the position information and the firstregistration; and controlling a display to display a portion of theultrasound image corresponding to the ultrasound imaging position and aportion of the medical image corresponding to the ultrasound imagingposition, based on at least one of the first registration and the secondregistration, wherein the performing of the first registration comprisesmatching the certain anatomical feature in the ultrasound image and thecertain anatomical feature in the medical image.
 16. The medical imagingmethod of claim 15, wherein the performing the first registrationcomprises: obtaining a first sample point of the certain anatomicalfeature extracted from the medical image of the object; obtaining asecond sample point of the certain anatomical feature extracted from theultrasound image of the object; and performing the first registrationbetween the medical image and the ultrasound image based on the firstsample point and the second sample point.
 17. The medical imaging methodof claim 15, wherein the ultrasound image is a three-dimensional (3D)ultrasound image, the medical image is a volumetric medical image whichincludes two-dimensional (2D) medical images obtained over a region ofthe object, the portion of the medical image includes one of the 2Dmedical images, and the display is controlled to display that one of the2D medical images which corresponds to the ultrasound imaging positionin the 3D ultrasound image, as indicated by the sensor.
 18. The medicalimaging method of claim 15, further comprising: obtaining a firstregistration information which provides a relationship between acoordinate system of the medical image and a coordinate system of theultrasound image based on the first registration; and calculating afirst transform between the coordinate system of the medical image andthe coordinate system of the ultrasound image, based on the firstregistration information, wherein the first transform is applied toperform a conversion between the coordinate system of the medical imageand the coordinate system of the ultrasound image.
 19. The medicalimaging method of claim 18, further comprising: performing a thirdregistration between the ultrasound image and the sensor, based on theposition information; and calculating a second transform which defines arelationship between a coordinate system of the sensor and thecoordinate system of the ultrasound image, wherein the performing thesecond registration comprises performing the second registration betweenthe medical image and the sensor based on a combination function of thefirst transform and the second transform.
 20. The medical imaging methodof claim 19, wherein the combination function is a product of the firsttransform and the second transform.